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Abstract 

Spike-train responses of single Hodgkin-Huxley (HH) and integrate-and-fire (IF) neu- 
rons with and without the refractory period, are calculated and compared. The HH and 
IF neurons are assumed to receive spike-train inputs with the constant interspike inter- 
vals (ISIs) and stochastic ISIs given by the Gamma distribution, through excitatory and 
inhibitory synaptic couplings: for both the couplings the HH neuron can fire while the 
IF neuron can only for the excitatory one. It is shown that the response to the constant- 
ISI inputs of the IF neuron strongly depends on the refractory period and the synaptic 
strength and that its response is rather different from that of the HH neuron. The vari- 
ability of HH and IF neurons depends not only on the jitter of the stochastic inputs but 
also on their mean and the synaptic strength. Even for the excitatory inputs, the type-I IF 
neuron may be a good substitute of the type-II HH neuron only in the limited parameter 
range 
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I. Inroduction 

A neuron generates the action potential, which propagates along the axon of a cell 
toward synapses exciting neurons in the next stage. There have been much theoretical 
studies on neurons which are responsible for encoding and decoding information carried 
by action potentials In the most popular neuron models like the Hopfield model 0], 
the output is described as a continuous variable which is slowly varying in time. The 
output is usually interpreted as short-time average of the rate of action potentials. This 
is based on the experimental observation that the mean firing rate depends on the applied 
stimulus to motor and sensory neurons. 

On the contrary, experimental evidences showing that the detailed timing and organi- 
zation of action potentials matter, have been reported in many biological systems; sonar 
processing of bats 0, sound localization of owls 0], electrosensation in electric fish |^, 
visual processing of cats [0, monkeys and human p. These suggest the importance 
of studying how neurons make computations based on the action potential timing with a 
resolution of the sub-millisecond range, receiving and emitting spike trains. 

Since Hodgkin and Huxley (HH)[0] first proposed the reliable model for squid giant 
axon, its property has been intensively investigated [|r^]-[|r7|. The HH model has been 
widely adopted for an investigation of biological systems with proper modifications of 



conductance channels |18|. Because the HH model is described by the non-linear differ- 
ential equations for four variables, its treatment is not easy, and then various types of 
simplified dynamical models have been proposed |jT9|. Among them, the simplest one is 
the integrate-and-fire (IF) model, which has been employed for a study on many kinds of 
subjects relevant to a single neuron as well as neural networks pO|. 



It has frequently claimed that the IF model captures the essentials of real neurons 



|2T|]-|2^. On the other hand, some studies have shown that the IF model is not realistic 
^-|2^. Actually the IF and HH models show important differences in their responses 
to applied, excitatory and inhibitory currents. For an excitatory (depolarizing) input 
dc current, I\, the IF neuron, which is classified as the type I, shows the self-excited 
oscillation with an arbitrary low firing frequency, /o [^. On the other hand, the HH 
neuron belonging to the type II has the discontinuous /o — I\ relation at the critical 
current, Jc, above which it shows the oscillation with a narrow range of /o (see Fig. 1(a)). 
For inhibitory (hyperpolarizing) input currents, the HH neuron may bring about a firing by 
the so-called rebound process against our intuition [^|, while the IF neuron cannot. This 
inhibitory rebound is realized in some non-linear dynamical neuron models like Fitzgh- 
Nagumo model. 

Quite recently, the present author has investigated the responses of the HH model 
to various types of spike-train inputs with constant, chaotic and stochastic interspike 
intervals (ISI): Ref.[^] is referred to as I in this paper. Our calculation in I shows that 
behavior of the variability to stochastic ISI inputs of the HH model is rather different 
from that obtained based on the IF model [^[^|. This is consistent with the recent 
calculation of Brown, Feng and Feerick, who show a stronger dependence of the 
variability on the level of inhibitory inputs in the IF neurons than in the HH neurons. 

It is the purpose of the present paper to elucidate the origin of the difference of the 
responses of IF and HH models to spike-train inputs, investigating whether the IF model 
may be an adequate substitute of the HH model. In order to clarify the point, we first 
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employ the spike-train inputs with simple constant ISIs, and then the stochastic ISI inputs 
to get some insight to the controversial variability problem initiated by Softky and Koch 

ii- 

Our paper is organized as follows: In the next sec. II, we describe a simple neuron 
model adopted for our numerical calculation. In sec. Ill, we investigate the response of our 
system to inputs with the constant ISI. Stochastic inputs with the Gamma distribution are 
treated in sec. IV. The final section VI is devoted to conclusion and discussion. Although 
some of the calculated results for the HH neuron have been published in I, we include 
them for the completeness of the present paper. 



2. Adopted model 

We adopt a simple system consisting of a neuron and a synapse; the former is described 
by the HH or IF model and the latter by the alpha function. We will investigate the 
response of our neuron when spike-train inputs are applied through the synapse. 

2.1 Hodgkin- Huxley model 

The HH model is described by the non-linear coupled differential equations for the 
four variables, V for the membrane potential, and m, h and n for the gating variables of 
Na and K channels, and it is given by [|10| 



CdV/dt = -g^.m^h{V - VNa) - 9Kn\V - Vk) - g^V - Vi) + (1) 

dm/dt = —{ajn + bm) rri + am, (2) 

dh/dt = -{ah + bh) h + ah, (3) 

dn/dt = —{an + &„) n + a„. (4) 

Here the reversal potentials of Na, K channels and leakage are V^a = 50 mV, Vk = —77 
mV and = —54.5 mV; the maximum values of corresponding conductivities are (^Na = 
120 mS/cm^, g^ = 36 mS/cm^ and gi, = 0.3 mS/cm^; the capacity of the membrane is 
C = 1 /iF/cm^; detailed expressions for am, bm et al. are presented in Refs.|T^ |31|] [0 . 



2.2 Integrate- and- Fire model 

We adopt the IF model which may include the absolute refractory period. The dy- 
namics of the membrane potential, V, and the phenomenologically introduced variable, 
p, is described by 

CdV/dt = -g{l+pa){V -V,-pVd) + {l-p)I, (5) 

dp/dt={-l/r,,)[p-e{p-w)], (6) 

with w = {Vt — V) / {Vt — V-r) , a = Tt./t^^ — 1 and Q{x) is the Heaviside step function. Here 
C is the capacitance and g {= C/rm) the conductance of the membrane; Xm is the life time 
of membrane potential and refractory period; and Vt are the reset and threshold 
potentials, respectively. In the limit of Tr ^ and Tp 0, our IF model given by Eqs.(5) 



and (6) is equivalent to the conventional IF model in which the reset condition is 



given by V{t^m) = ^ V{t'^m) = at tom , the firing time of \^ (m = 0, 1, 
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The variable p is zero in the active period and it is one in the refractory period. The 
rapid transient between the two states of p is given by Eq.(6) with a small Tp. In the 
active period (p = 0), Eq.(5) becomes 

dV/dt = -{l/r^){V - + //C, (7) 



which is just the same as that given in the conventional leaky IF model |]20[. On the other 
hand, in the refractory period (p = 1) Eq.(5) becomes 

dV/dt = -(l/r,)(V -V,- Ki), (8) 

where Tr = t^/{1 + a). Equation (8) shows that the reset of V is accomplished by the 
deriving potential, + with the time constant t,. which is smaller than t^. 

The advantages of the present model given by Eqs. (5) and (6) are that it includes the 
refractory period and that it automatically resets the membrane potential. The model 



similar to ours was previously proposed by Horn and Opher|]5S[. The physical meaning of 
their model is, however, not transparent and their a term, which corresponds to our Vd 
term in Eq.(6), improperly persists in the active period. It is crucial to reject this term in 
the active period (Eq.(7)) for a discussion of the input-output response of the IF neuron. 

The bold curve in Fig.l expresses the /o — /j relation of the HH neuron, showing the 
discontinuous transition at 1^ = 6.3/iA/cm^. It is rather difficult to choose the parameters 
of the IF model so as to reproduce the fo — Ii relation of the HH neuron because the IF and 
HH neurons belong to the different type neurons, although the exact choice of parameter 
values is not crucial. After several tries |^6[, we have determined to employ the following 



parameters for our numerical calculations: = —75 mV, Vt = —55 mV, Vd = —10 mV, 



Tm = 20 msec Tp = 0.02 msec and C = 4yuF/cm^. As for the refractory period, t,., 
we adopt the two choices; neurons with Tr = 0.1 and 2.0 msec are referred to as the IFO 
and IFl neurons, respectively. 

After a simple calculation, we get the fo dependence on Ji of IF neurons, given by 

l//o = To = Ht^) + ln[ ^"~^,^^~^^^ ] + 0(rp), (9) 

— Vd 

where the critical current, Jc, is given by 

/c = (C/rJ(V;-K), (10) 

leading to Jc = 4 /iA/cm^. Dashed and solid curves in Fig.l express fo for IFO and IFl 
neurons, respectively. 

Figure 1(b) shows the examples of the self-excited oscillations of the membrane po- 
tentials of the HH, IFO and IFl neurons. For IF neurons we plot V + cp{l — p) instead 
of V: the c term (c = 350) yields the spiky contribution. 

2.3 Synaptic inputs 

We consider the delta- function- type spike-train input expressed by 

f^i(i) = KE^(^-U, (11) 
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where is the magnitude of the action potential, and the firing time ti„ for arbitrary 
is assumed to be recurrently defined by 



n 



til = 0, (13) 

ISI of input spike, Ti„, being generally a function of a given time ti„. 

It has been reported that biological synapses exhibit temporal dynamics during neu- 



ronal computations [^] [37|. We, however, treat the synapse as a static unit for a simplic- 
ity of our calculation. The spike train given by Eq.(ll) is injected through the synapse, 
yielding the postsynaptic current Ji given by 

/i(t) = ^ a(t-ti„), (14) 

n 

where Asyn = Qsyn (Va — Vsjn), f/syn and V^yn are the synaptic conductance and reversal 
potential, respectively, and the alpha function, a(t), is defined by 



«(t) = (t/r,) e-*/^= e(t), (15) 

Tg being the time constant relevant to the synapse conduction. The positive and negative 
Agyn stand for the excitatory and inhibitory couplings, respectively. When ^ r^, 
Eqs.(14) and (15) yield pulse currents with the maximum value of Jj™^^ = 0.368 Asyn at 
t = Tin + Ts and with the half-width of 2.45 Tg. On the contrary, when ^ Tg the 
temporal summation of input currents is realized because an input pulse comes before the 
current induced by its preceding pulse is not attenuated. We assume = 2 msec and 
treat A^y^ as an adjustable parameter. 

When the membrane potential V oscillates, it yields the spike-train output, which may 
be expressed by 

Uo{t)=Va E 5(t-to™), (16) 
m 

in a way similar to Eq.(ll), where tom is defined as the time when the membrane potential 
V{t) crosses = mV from below. The output ISI is defined by 

Tom tom+1 torn- (-^'^) 

Differential equations given by Eqs.(l)-(4) for the HH neuron (or Eqs.(5) and (6) 
for the IF neuron) including the external current given by Eqs.(ll)-(15) are solved by 
the forth-order Runge-Kutta method with the integration time step of 0.01 msec. The 
calculations are performed for 2 sec in the constant-ISI case and for 20 sec in the stochastic 
ISI case. If ISI of spike-train input or output is about 10 msec in the latter case, the size 
of its sample is about 2000. Although this figure is not sufficiently large for statistics of 
ISI data, we hope an essential ingredient will be clarified in our numerical investigation. 
Hereafter time (Ti„ etc.), voltage {V etc.) and current (Ji, A^yn etc.) are expressed in 
units of msec, mV and /iA/cm^, respectively unless otherwise specified. 



3. Constant-ISI inputs 
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3. 1 HH neurons 

Let us first consider the HH neuron which receives spike-train inputs given by Eqs.(ll)- 
(15) with the positive Agyn and constant Ti„. The sohd curve in Fig. 2(c) expresses the 
time course of the membrane potential, V, for the excitatory postsynaptic current, /i, with 
^syn = 40 shown by the sohd curve in Fig.2(b), which is induced by an apphed spike-train 
input, Ui, with Ti„ = 20 msec as depicted in Fig.2(a) (for time courses of potentials for 
inputs with Ti„ = 10 msec, see Fig. 5 of I). For this spike-train input of Ti„ = 20 msec, we 
get a regular spike-train output with Tom = 20 msec. This is in contrast with the case of 
spike-train input of Ti„ = 10 msec, in which output spike-train is phase locked with the 
ratio of 4 : 3, oscillating with a long cycle of 40.00 msec (=11.25 + 12.36 + 16.39) = 4Ti„, 
where 11.25, 12.36 and 16.39 are the values of output ISIs. 

Solid and dashed curves in Fig. 3 (a) show the average (/io) and root-mean-square (RMS, 
(To) values of {Tom} for input ISI of Ti„ = 10 msec when A^yn is changed . Filled circles 
denote the distribution of {Tom} for a given Agyn- For example, in the case of Asjn = 40, 
we get Tom = 11.25, 12.36 and 16.39 msec as mentioned above. For Agyn > 56, ISI of the 
output is the same as that of input, then the coefficient defined by = /io//^i is unity. 
When Asyn is decreased, we get k larger than unity, and /c = 2 for 8 < A^y^ < 28. Below 
^syn = 8, we have no output spikes. 

Solid and dashed curves in Fig. 3(b) express /io and ctq when the average of input ISI 
(/ii) is changed with Agyn = 40. We obtain k = 1 for /i; ~ 12 msec, = 2 for 6 ~ yUj ~ 8 
msec, and k = 3 for fii = 4 msec: otherwise k is non-integer. We obtain no output ISIs of 
Tom < 10 msec for Ti„ < 12 msec, which is due to the low-pass filter character of the HH 
neuron ||31[| . 

Next we discuss the inhibitory input case. The dash curve in Fig. 2(c) expresses the 
time course of membrane potential with Tom = 20 msec when we apply the inhibitory 
postsynaptic current with Asyn = —40 and Ti„ = 20 msec, which works to hyperpolarize 
the membrane potential, V. When the postsynaptic current given by the alpha function 
(Eq.l5) decreases and vanishes, V changes to restore to the rest level and it crosses the 
threshold to yield an action potential with a delay of about 15 msec. Thus the HH neuron 
can fire even for the inhibitory input by the rebound process This process, however, 
requires an appreciable periods: when the input ISI is less than about 15 msec, the HH 
neuron cannot fire because the next inhibitory spike is applied before the hyperpolarized 
membrane potential crosses the threshold level. 

By changing Agyn and /ii, we perform similar calculations, whose result is summarized 
in Fig.4(a). It expresses the phase diagram showing the region where the integer k is 
obtained in the /ii — Agyn space. Note that between the integer fc's, we have non-integer 
solutions; for example, in the case of /Zi = 10 msec and Asyn = 40, we get k = 1.33 between 
k = 1 and 2. We notice for positive Agyn that there is the wide k = 1 region and that the 
regions of larger k appear at the left side of the k = 1 region. For negative Agyn, the HH 
neuron fires for input with /i > 15 msec, for which we get the states of A; = 1 and k = 2. 

3.2 IF neurons 

Now we consider the IF neuron to which the constant-ISI input is applied. The solid 
(dashed) curve in Fig.2(d) shows the time course of the membrane potentials of the IFl 
neuron for excitatory (inhibitory) inputs. When the excitatory input with Ti„ = 20 msec 
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is applied, its membrane potential is depolarized to cross the threshold, and it emits 
the spike-train output with Tq^ = 20 msec. When the inhibitory input is applied, its 
membrane potential is hyperpolarized. After postsynaptic current vanishes, V changes to 
restore to the rest level but cannot cross the threshold level. The behavior of V of the IFl 
neuron is in strong contrast with that of the HH neuron shown by the dashed curve in 
Fig.2(c). Thus the IFl (and IFO) neuron is not excitable for inhibitory spike-train inputs. 

Figure 5 shows the time courses of Ui{t) and V{t) of the IFl neuron for Ti^ — 10 msec 
with Agyn = 40, 64 and 120. For Agyn — 64, IFl neurons regularly emit spike trains of 
Tom = 10 msec with a delay of about 3 msec. For Agy^ = 40, the output ISI becomes 
larger than 10 msec, because it can emit output pulse after integrating small input signals 
by virtue of the type-I neuron although a single input pulse is insufficient to trigger the 
output pulse. On the contrary, for a large A^y^ = 120, a single input pulse make the IFl 
to emit irregularly multiple pulses with Tom smaller than Tin- The time course of the 
membrane potential of the IFO neuron is ostensibly similar to that of the IFl neuron (not 
shown) . 

The ^syn dependence of the responses of the IFO and IFl neurons to inputs with 
l^i = 10 msec is shown in Fig.6(a) and 7(a), respectively. Solid (dashed) curves express 
Ho (co) and filled circles denote {Tom} for a given Agya- For the IFO neuron, we get the 
integer values of k — 1 for 51 ~ Agy^ ~ 58, A; = 2 for Agy^ — 32, k — 3 for Asyn — 26 
and k = 4 for A^y^ = 22. On the contrary, for the IFl neuron we get the integer values 
of A; = 1 for 54 ~ Asyn ~ 83, A; = 2 for 32 ~ Asyn ~ 37, A; = 3 for A^yn = 26 and A; = 4 
for Asyn = 22. When we compare these results with that of HH neurons (Figs. 3(a) and 

(b) ), we notice that A; ~ 3 is realized for small A^y^ and that A; < 1 exists at large Agyn. 
Figures 6(b) and 7(b) show the /^i dependence of //q and tJo of the IFO and IFl neurons, 

respectively, with Agyn — 64. We get //o — Hi {k — 1) for //i > 15 (//i > 6) msec in the IFO 
(IFl) neuron. When Agyn is reduced, we obtain the states with larger k. For example, for 
Asyn = 40, we get = 2 with 16 < < 22 (12 < Hi < 22) msec in the IFO (IFl) neuron; 
the A; = 1 state is not available (not shown). 

From calculations by changing Asyn and /Xi, we obtain the phase diagrams of the integer 
k values in the for the IFO and IFl neuron, which are shown in Figs.4(b) and 

(c) , respectively. Because IFO and IFl neurons cannot fire for inhibitory inputs, results 
only for the positive Asyn are shown. A comparison of them with the corresponding phase 
diagram of the HH neuron (Fig. 4 (a)) shows 

(1) the A; = 1 region, particularly of the IFO neuron, is greatly reduced, 

(2) although the k > 2 region in the HH neuron appears at the left side of the A; = 1 
phase, such phases in IFO and IFl neurons appear below the k = 1 phase, 

(3) the phase with A' < 1 appears for IFO and IFl neurons although it does not exist for 
the HH neuron in the parameter range shown in the Fig.4(a), and 

(4) an inclusion of the refractory period in the IFl neuron widens the A; = 1 region. 

4. Stochastic-ISI inputs 

The ISIs of spike-train input. Tin, in Eq.(12) are assumed to be independent random 
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(18) 



for which we get /li — r/s, and (7i = y/r/s, F (r) being the gamma function. For a later 
purpose, we define the dimensionless variabihty given by 



for input (A = i) and output ISIs (A = o), from which we get Cvi = l/v^ foi' the Gamma- 
distribution inputs. For r = 1 in Eq.(18) we recover the exponential distribution (cvi — 1) 
and a Poisson distribution for the number of spikes in a given time interval. In the limits 
of r — s> oo and s — cx) with keeping fii = r/s fixed, Eq.(18) reduces to P{T) — S{T — /li), 
the constant ISl with /ii = T and Cvi = 0. 

The spike-train input created by the Gamma-distribution generator is applied to our 
neural system. Calculations are performed by changing Asyn or /li by keeping the value 
of Cvi fixed. Note that because the number of our sample of input ISI is not sufficiently 
large, the obtained Cyi fluctuates around the intended values. 

4-1 HH neurons 

The histogram in Fig. 8(a) shows the distribution of excitatory input ISIs (Agyn = 40) 
with Cvi = 0.40 and jJLi — 10 msec, which leads to output ISIs with the distribution shown 
in Fig.8(b) (cvo = 0.25 and = 14.84 msec). Output histogram in Fig.8(b) has no 
distributions at Tom < 10 msec, which arises from the low-pass filter character of the HH 
neuron as shown in Fig. 3(b). 

Figures 9(a) and (b) show the results of //q, (Jq a-nd Cyo as a function of A^yn for inputs 
with Cvi = 0.4 and 1.0, respectively (//i = 10 msec). Comparing them with the results 
for the constant ISI (cyj = 0) shown in Fig.3(a), we note that for larger c^i, values of 
/io > 20 msec are much accumulated in the region of smaller Agyn region. This is due 
to the integration character of neurons for weak inputs and it is also realized in the IF 
neuron (Figs. 10(a) and 11(a)). 

The dependence of output ISIs on ji^ for Cvi = 0.4 and 1.0 are shown in Fig.9(c) and 
(d), respectively. We note that the structure at < 10 seen in the case of Cvo = 
(Fig.3(b), disappears because of the randomness in inputs. Although fi^ ~ for large 
lii, we get > 10 msec for small fii, since the HH neuron plays as the low-pass filter. 
For Cvi = 1.0, we get ~ (/Xi -|- 10) whereas cTq ~ /^i, which yield an increase in Cvo as 
increasing /li. 

Figures 9(e) and (f) express the /ij dependence of output ISIs for inhibitory inputs 
with Cvi = 0.4 and 1.0, respectively. Their comparisons with the result for Cvi = 0.0 (Fig. 3 
(c)) show that the HH neuron may fire for stochastic inputs with fii < 15 msec. It is 
interesting to note from a comparison of Fig. 9(c) (Fig.9(d)) with Fig.9(e) (Fig. 9(f)) that 
the /Ii dependence for inhibitory inputs is quantitatively similar to that for excitatory 
inputs. 

4-2 IF neurons 

Histograms in Figs. 8(c) and (d) show the distributions of output ISIs of the IFO neuron 
and those of the IFl neuron, respectively, for excitatory inputs of Cvi = 0.4 and /ii = 10 



CvA ^ CTx/lJ'X 



(A = i and o). 



(19) 
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msec whose distribution is plotted in Fig. 8(a). Solid and dashed histograms in Fig.8(c) 
and (d) express the results of Agyn = 40 and 64, respectively. The input histogram has 
a peak at about Tin = 10 msec, as expected. On the contrary, output ISls of IFO and 
IFl neurons for Asy^ = 64 have peaks at lower values of ISI. This is understood from 
Figs. 6(a) and 7(a) which show that A; < 1 for fi[ < 16 (12) msec in the IFO (IFl) neuron. 
Furthermore, for the IFl neuron, the refractory period makes the neuron to emit no 
outputs for Tom < 4 msec. 

Figures 10(a) and (b) (11(a) and (b)) show output results for excitatory inputs with 
Cvi = 0.4 and 1.0, respectively, of the IFO (IFl) neuron as a function of Agy^. As Cvi is 
increased, large /Iq values are accumulated in the region with a small Asyn- Although this 
phenomenon is realized also in the HH neuron (Fig.9(a) and (b)), the fio value is 100 
msec at most in the HH neuron while the maximum /Xq value exceeds 1000 msec in the IF 
neurons. The difference in /Iq and cxo between the IFO and IFl neurons become smaller 
for larger Cvi, although Cvo of the IFO neuron is always larger than that of the IFl neuron. 

The /ii dependence of the output ISIs for excitatory inputs with Cvi = 0.4 and 1.0 of 
the IFO (IFl) neuron are shown in Fig.lO(b) and (d) (11(b) and (d)), respectively. With 
increasing /Xj, both /ij and ai increase, but Cvo tends to saturate. For larger q, we get 
larg Gr Co 5 3bS expected. It is interesting to compare the results of the IF neurons with 
those of the HH neuron. For inputs with Cj = 0.4 and /ii = 10 — 30 msec, we obtain 
Cvo = 0.25 - 0.40, 0.47-0.77 and 0.40-0.45 for the HH, IFO and IFl neurons, respectively. 
Similarly, for inputs with q = 1.0 and /ii = 10 — 30 msec, we get Cvo = 0.56 — 0.86, 
1.01-1.13 and 0.85-0.96 for the HH, IFO and IFl neurons, respectively. We note that the 
relation: 

cr < ci^^ < cr - c.i, (20) 

holds in our calculation, related discussion being given in the next section. 



5. Conclusion and discussion 



Since Softky and Koch ||33[ reported a large Cvo (~ 0.5 — 1.0) in cortical neurons in 
visual VI and MT of monkeys, it has been controversial how to understand the large 
variability of neurons to stochastic inputs There have been much discussions 



on this subject using the IF model. Some theoretical studies show that IF models lead to 
small Cvo because an integration of a large number of random inputs works to reduce the 
variability On the other hand, other studies have shown that IF neurons may 

yield an appreciable value of Cvo |^ • 

Figure 12 shows the Cvi — Cvo plot of ISI data having reported for HH (Figs. 3(b), 9(c) 
and 9(d)), IFO (Figs.6(b), 10(c) and 10(d)) and IFl neurons (Figs.7(b), 11(c) and 11(d)) 
and of new results calculated for Cvi = 0.7. Open squares, filled triangles and circles 
denote the results of the HH, IFO and IFl neurons, respectively. We note the variability 
of output ISIs may be large and nearly the same as that of input ISIs in HH and IF 
neurons. Scattered values of Cvo for a given Cvi mean that Cvo depends not only on Cvi but 
also strongly on /ii (and Agyj^). For example, when the /ii is varied from 2 to 30 msec in 
the IFO neuron, we get c^ = 0.20 - 0.48 for c^i = 0.4, Cvo = 0.40 - 0.84 for Cvi = 0.7 and 
Cvo = 0.64 — 1.18 for Cvi = 1.0. Our calculation reconciles the dispute among the earlier 
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calculations yielding small |2^[Q and large Cvo|^|2S||lT| based on the IF model; the 
difference in the calculated Cvo may be due to the difference in the parameters adopted in 
their calculations. 

It should be noted in Figs. 9, 10 and 11 that the variability of output ISIs generally 
increases with increasing input ISIs, which is in agreement with the biological data (Fig. 3 
of Ref.|Q). When the input ISI is small compared to the characteristic integration time, 
the neuron acts as an integrator yielding a small Cvo, while when the reserve is true, the 
neuron play a role of the coincidence detector with a fairly large Cvo- 

To summarize, we have performed numerical calculations of the responses of the HH 
and IF neurons to inputs with the constant and stochastic ISIs, to make a comparison 
between them, after we had chosen the parameters such as for the IF model to mimic 
the time dependence of the membrane potentials of the HH neuron (Fig. 1(a)). Our 
calculations have shown the followings: 

(i) the HH neuron can fire for both excitatory and inhibitory inputs whereas the IF neuron 
only for excitatory inputs, 

(ii) responses of the HH and IF neurons show the complicated behavior to spike-train 
inputs even with the simple constant ISIs, for which /io (or k = Ho/ni) is generally 
functions of /ii, Cvi and Agyn, 

(iii) the response to constant ISIs of the IF neurons is rather different from that of the 
HH neuron: the k = 1 region of the IF neuron, in which neurons properly respond to 
inputs, is much narrower than that of the HH neuron (Fig. 4), 

(iv) the variability of IF and HH neurons to stochastic ISIs shown in Fig. 12, may be large 
and follow the relation given by Eq.(20), 

(v) the Asyn dependence of the variability of the IF neuron is stronger than that of the 
HH neuron, and 

(vi) an inclusion of the refractory period in IFl depresses the Cvo values and improves to 
some extent its response, widening the k = 1 region. 

The difference in the item (i) arises from the fact that the HH model has the rebound 
process against the hyperpolarized membrane potential while the IF model does not. In 
order to make the item (ii) more concrete, we show in Figs. 13(a) and (b), the f\ — fo plots 
of the HH and IFO neurons, respectively, where f\ (= l/fJ^x) is the mean frequency of 
input (A = i) and output ISIs (A = o). Saturating functions with the maximum frequency 
of /™^^ ~ 100 Hz given by the /i — fo plots of the HH neuron (Fig. 13 (a)), are similar to 
the sigmoidal function of g{x) = (1 + tanha;)/2 adopted in the formal rate-coding models 
0. However, monotone increasing functions given by those of the IFO neuron (Fig. 13(b)) 
are quite different from g{x). The /i — fo plots of the IFl neuron express also increasing 
functions but saturate with the large maximum frequency of f^'^^ ~ 500 Hz (not shown). 
The main origins yielding the differences cited in the items (iii)-(vi) are (1) the difference 
of the fo — Ii relation of the type-I IF neuron from that of the type-II HH neuron and (2) 
the difference in the refractory period. The continuous fo — Ii relation in the IF neuron 
yields a large k {> 3) after integrating small inputs for a small Agyn- The IF models 
with the vanishing or small refractory periods emits the output ISI with A; < 1. Then the 
k = 1 region in the phase diagram of the IF neuron is much reduced compared with that 
of the HH neuron (Fig.4). The item (iv) is consistent with earlier calculations using the IF 
|2^|2^[^ and HH models |2^. Brown, Feng and Feedick show that the variability 



11 



of the IF neuron has a stronger dependence on the number of synaptic inputs, A^s, than 
that of the HH neuron. Since their A^s is expected to correspond to our Asyn in a crude 



sense, the item (v) is consistent with their result. They also claim that an inclusion 
of the absolute refractory period in the IF neuron decreases the variability and that it 
increases the disparity between the results for the IF and HH neurons. The former agrees 
with our item (vi) but the latter does not. Gutkin and Ermentrout predict based on 
the Morris-Lecar (ML) model that the variability of type-I ML neuron is larger than that 
of type-II ML one, which is supported by our calculations (Eq.20)). However, their claim 
that the type-II neuron yield only a small Cvo, cannot be applied to the HH model. 

The HH neuron responds to static and spike inputs differently from the IF neuron. The 
HH neuron shows the complex behavior not shared with an IF neuron. The chaotic oscil- 
lation induced by an applied sinusoidal currents [00] and the firing by the inhibitory 
rebound are never realized in the single IF neuron. Our calculations have shown that 
even for the excitatory inputs the IF model may be a good substitute of the HH model 
only within the limited parameters. There are many experimental and theoretical evi- 
dences showing that reciprocally inhibitory neurons play important roles in real systems 
such as Hippocampus and thalamus for which the IF model cannot be used. We 
should mind advantage and disadvantage of the IF neuron in modeling biological neural 
systems. 
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Figure Captions 

Fig.l (a) The fo—Ii plot of the HH (bold solid curve), IFO (dashed curve) and IFl neurons 
(solid curve), and (b) the self-excited oscillations of the HH, IFO and IFl neurons, scales 
for IFO and IFl neurons being shifted by 200 and 400 mV, respectively. 

Fig.2 The time courses of (a) the input spike train, Ui, with Tin — 20 msec, (b) the 
postsynaptic currents, /j, and (c) the membrane potential, V, of the HH (Agyn = ±40 
/iA/cm^) and (d) V of IFl neurons (Agyn = ±64 /xA/cm^); solid and dashed curves 
denote excitatory and inhibitory cases, respectively, and scales for Ui and /i are arbitrary, 
(see text). 

Fig.3 Responses of the HH neuron to constant ISIs. (a) mean (/Xq, solid curve), RMS 
((To, dashed curve) and the distribution (filled circles) of output ISIs as a function of Asyn 
for = 10 msec (the arrow denotes the input ISI); (b) those as a function of fi^ with 
^syn = 40 /xA/cm^ and (c) with Agyn = —40 iiK/cw? (dotted lines express k = 1, 2 and 
3). 

Fig.4 The phase diagrams of the (a) HH, (b) IFO and (c) IFl neurons in the fii—Agyn space, 
showing the states with integer k (= /Xo/yLii) for constant-ISI inputs given by Eqs.(14) and 
(15), crosses denoting no outputs. The results of IFO and IFl are shown only for positive 
Agyn because they cannot fire for inhibitory inputs. 

Fig.5 The time courses of (a) constant-ISI input and (b) the membrane potential of IFl 
neuron with Agyn = 40, (c) Agyn = 64 and (d) Asyn = 120 /xA/cm^. The scale of (a) is 
arbitrary and those of (c) and (d) are shifted by 200 and 400 mV, respectively. 

Fig.6 Responses of the IFO neuron to constant ISIs. (a) mean (//q, solid curve), RMS (uo, 
dashed curve) and the distribution (filled circles) of output ISIs as a function of Asyn for 
^i=10 msec (the arrow denotes the input ISI); (b) those as a function of /li (dotted lines 
express k = /Xo//Xi = 1). 

Fig.7 Responses of the IFl neuron to constant ISIs, same as in Fig.6. 

Fig. 8 Histograms of (a) the stochastic input ( /ii = 10 msec, Cvi = 0.4), (b) output ISIs of 
the HH neuron, (c) of the IFO neuron and (d) of the IFl neuron, solid (dashed) histograms 
being for Agy^ = 40 (64) /xA/cm^. 

Fig.9 Responses of the HH neuron to stochastic ISIs; (a) /Iq (solid curve), cTq (dashed 
curve) and Cvo (thin solid curve) as a function of Agyn for inputs of yUi = 10 msec with 
Cvi = 0.4 and (b) with Cyi = 1.0 (the arrow denotes /Xj); (c) those as a function of /xi with 
^syn = 40 /xA/cm^ for inputs of Cvi = 0.4 and (d) of Cyi = 1.0; (e) those as a function of 
/Xi with Asyn — —40 /xA/cm^ for inputs of Cyi = 0.4 and (f) of Cyi = 1.0. 

Fig. 10 Responses of the IFO neuron to stochastic ISIs; (a) /io ( solid curve), ag (dashed 



15 



curve) and Cvo (thin solid curve) as a function of Agyn for inputs of /^i = 10 msec with 
Cvi = 0.4 and (b) with Cvi = 1.0 (the arrow denoting /xi); (c) those as a function of /li with 
^syn = 64 nA/cm^ for inputs of Cvi = 0.4 and of Cvi = 1.0. 

Fig. 11 Responses of the IFO neuron to stochastic ISIs, same as Fig. 10. 

Fig. 12 Cvo against Cyi of HH (open squares), IFO (triangles) and IFl neurons (circles) (see 
text). 

Fig. 13 The /j — fo plot of (a) the HH and (b) IFO neurons, dashed curves denoting the 
extrapolation. 



This figure "Figl-4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9907387v4 



This figure "Fig5-8.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9907387v4 



This figure "Fig9-10.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9907387v4 



This figure "Figll-13.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9907387v4 



